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Using a non-perturbative functional renormalization group approach we calculate the renormalized 
quasi-particle velocity v{k) and the static dielectric function e(fc) of suspended graphene as functions 
of an external momentum k. Our numerical result for v{k) can be fitted by v{k)/vF = A-\-B ln(Ao/fc), 
where vf is the bare Fermi velocity, Aq is an ultraviolet cutoff, and A — 1.37, B = 0.51 for the 
physically relevant value (e^ jvF = 2.2) of the coupling constant. In contrast to calculations based 
on the static random-phase approximation, we find that tik) approaches unity for A: —> 0. Our result 
for vik) agrees very well with a recent measurement by Elias et al. [Nat. Phys. 7, 701 (2011)]. 
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At low energies the physical properties of graphene 
are dominated by the Dirac points where the energy 
dispersion vanishes linearly. In this regime many-body 
effects become important and can be measured experi¬ 
mentally [1]. In view of the great interest in graphene 
both for fundamental research and applied physics, it is 
important to gain a thorough understanding of correla¬ 
tion effects. Of particular interest is the renormalization 
of the Fermi velocity at the Dirac points by long-range 
Coulomb interactions, which has been observed exper¬ 
imentally in suspended graphene using cyclotron reso¬ 
nance [2], in ARPES measurements of quasi-freestanding 
graphene on SiC [3], and in graphene on hexagonal boron 
nitride (hBN) [4]. Early one-loop renormalization group 
(RG) calculations [5] predicted a logarithmic enhance¬ 
ment of the renormalized Fermi velocity, 

vk/vf = l + {a/A)hi{KQ/K)FO{a^), (I) 

where A is the infrared cutoff introduced in the RG proce¬ 
dure, Aq is an ultraviolet cutoff of the order of the inverse 
lattice spacing, vp = 10®m/s is the bare Fermi velocity, 
and a = jvp is the relevant dimensionless coupling 
constant. Because for graphene suspended in vacuum 
a « 2.2 is rather large, perturbative RG calculations are 
not expected to be quantitatively accurate. 

In this work, we use a functional renormalization group 
(FRG) approach [6, 7] to derive non-perturbative RG 
flow equations for the cutoff- and momentum-dependent 
velocity v\{k) and the static dielectric function eA{q) of 
suspended graphene. Since we are interested in the RG 
flow of momentum-dependent functions, the field theo¬ 
retical RG is not sufficient, because with this method 
one can only keep track of a finite set of coupling con¬ 
stants. We show here that this problem can be solved 
within the FRG formalism [6, 7]; specifically, we de¬ 
rive two coupled integro-differential equations for the 
cutoff-dependent functions ua(A:) and eA{q) which are 
non-perturbative in a and self-consistently describe the 
interplay between self-energy and screening effects. 


Our starting point is the following effective Hamilto¬ 
nian describing the low-energy physics of graphene, 

= X + A [ fqP-qPq, (2) 

P Jk Jq 

where p = ± labels the two Dirac points of the un¬ 
derlying tight-binding model on a honeycomb lattice, 
Vp = pvp is the bare Fermi velocity at Dirac point p, 
and tpp^k) are two-component fermionic held operators 
whose components are associated with the two sublat¬ 
tices of the honeycomb lattice. The two-component vec¬ 
tor cr = (ct“, (t^) contains Pauli matrices acting on sublat¬ 
tice space, and two-dimensional momentum integrations 
are denoted by /^^ = / The interaction in Eq. (2) is 

specihed in terms of the Eourier transform fq = 27re^/|<7| 
of the Coulomb interaction and the Eourier components 
of the density operators, pq — J2p Ik + Tor 

simplicity, we consider a given spin projection and sup¬ 
press the spin label. We shall insert the spin-degeneracy 
= 2 of the electrons in Eq. (10) below. 

To derive FRG how equations, we introduce a cutoff 
A which inhibits the propagation of electrons with mo¬ 
menta I A; I < A. For our purpose it is sufficient to work 
with a sharp momentum cutoff. The regularized free 
propagator is then Gp a(A') = 0(fc — A) [iuj — Vpcr ■ k]~^, 
where the label K = (fc, iuj) represents momentum k and 
fermionic Matsubara frequency iuj. At some large initial 
cutoff Aq of the order of the inverse lattice spacing the 
regularized Euclidean action of our system is 

p Jk 

+ ^^p{-QmQ)], ( 3 ) 

where 'ipp{K) is a two-component Grassmann held and 
we have used a Hubbard-Stratonovich transformation to 
represent the Coulomb interaction in terms of a scalar 
held 4>{Q) which couples to the Eourier components 
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P{Q) = Ep/k + Q) of the density. The 

integration symbols are /^ = /^ / ^ and /^ = /^ / g. 
Here Q = {q,ioj), where iCo is a bosonic Matsubara fre¬ 
quency. It is now straightforward to write down a for¬ 
mally exact flow equation for the generating functional 
rA['i /'5 </*] of the irreducible vertices describing their change 
as we reduce the infrared cutoff Ao —> A. By construc¬ 
tion for A —)■ 0 the flowing vertices reduce to the exact 
irreducible vertices of our original Hamiltonian (2). As 
usual, we obtain an approximate solution of this func¬ 
tional flow equation by working with a truncated form 
of rA[f/') For our purpose, the following truncation is 
sufficient, 

p ^ 

+ l [ H-Q)F^\Q)HQ) 

J Q 

+ E/ [ rlAiK,Q)^p;iK + Q)rpiK)^{Q), (4) 

p,s JkJq 

where in the last term s = A, H is the sublattice la¬ 
bel and ipp{K), 'ipp{K) are the sublattice components of 
'tpp{K) and the adjoint spinor It is convenient to 

express the renormalized fermionic and bosonic propaga¬ 


tors in terms of the corresponding self-energies as usual, 
G-iW = [G%(^)]-'-Sp,A(iF), and F^\Q) = f-^ + 
Ha(< 5). Note that the fermionic self-energy T,p^A{K) is 
a matrix in the sublattice labels. Our truncation (4) re¬ 
tains only those vertices which are already present in the 
bare action (3). Although higher order vertices with more 
than three external legs are generated by the RG proce¬ 
dure, they are irrelevant at the Gaussian fixed point. In 
fact, a simple scaling analysis keeping the Gaussian part 
of rA[f/', 4^\ invariant shows that the renormalized vertices 
r„ with n external legs scale as A"“^, implying that all 
vertices with n > 3 external legs are irrelevant. There 
are hve marginal vertices with three external legs, cor¬ 
responding to the field combinations 

(j)^ and (jxfxj). While in our cutoff scheme 

the purely bosonic (/)0^-vertex does not couple to the 
flow of the self-energies F,p{K) and n((5), the sublattice¬ 
changing vertices of the type and (p are ne¬ 

glected in Eq. (4). We nevertheless believe that our trun¬ 
cation (4) is accurate because the sublattice-preserving 
vertices of the type {p^ip^cp and 'tp^'ip^(p are already h- 
nite in the bare action (3). 

Within our truncation and cutoff scheme, the self¬ 
energies and three-legged vertices appearing in Eq. (4) 
satisfy the following system of ERG flow equations. 


dA^fiK) 

dAli{Q) 

dAr;{K,Q) 


F(g)G- {K - Q)TUK - Q, Q)r! {K,-Q), 


EE/J Gf{K)G^p^K -Q) + Gf{K)G‘p%K - Q)\r;{K, -Q)r‘p{K -Q,Q), 

ss' 

I F{Q') I Gl^' {K + Q- q')g;%k - g') + Gf {K + Q- g')Gr (iF - g') 


’Q' 


xr;{K + Q-Q', -gOr;' {k - g', g). 


(5a) 

(5b) 

(5c) 


where for simplicity we have omitted the cutoff label A 
and the single-scale propagator is 

Gp,a{F) = -S{k - A) [zw - VpCT ■ k - Ep^AiK)]~^. (6) 

The derivation of the above ERG equations starting from 
the general Wetterich equation [8] for the theory defined 
in Eq. (3) is similar to the derivation of the vertex expan¬ 
sion for mixed Bose-Fermi models discussed in Refs. [6 
and 9]. A graphical representation of the flow equa¬ 
tions (5a-5c) is shown in Fig. 1. Since we are interested 
in the low-energy behavior of the Green function, we ex¬ 
pand the self-energy to linear order in the frequency, 

^p,a{K) = pVA{k)a- • fe + (1 - Zpp^)iio + 0(cc2). (7) 

Here Za is the wave-function renormalization factor. 
Note that we do not expand the momentum dependence 


of the velocity correction VA^k). For small frequencies 
our scale-dependent fermionic propagator is given by 


Gp,a{K) 


-e(fc- a)Za 


iuj + pvA{k)(T ■ k 


( 8 ) 


where ^A(fc) = VA{k)k is the energy dispersion at cutoff 
scale A and VA{k) = Za[vf F VA{k)] is the correspond¬ 
ing velocity. Moreover, we retain only the marginal part 
E® a(0;0) = *7A of the three-point vertices. The flow 
equation (5a) for the self-energy then reduces to 


Oa^pAK) 



Sjq-A fk-q 
eAik — q,i^ — fw) 


iu) +pvA{q)o' ■ q 


(9) 


Here we have introduced the cutoff-dependent dielectric 
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+ ^ '^''O*^ 


^ A = 

B • B 


FIG. 1. (Color online) Diagrammatic repesentation of the 
FRG flow equations (5a-5c) for the fermionic self-energy (first 
two lines), for the polarization (third and fourth line), and 
the three-legged vertices (last two lines). Solid arrows repre¬ 
sent the exact cutoff-dependent fermionic propagators, wavy 
lines represent the corresponding bosonic propagators, and 
single-scale propagators have an additional slash. The dot 
over the vertices on the left-hand side represents the deriva¬ 
tive with respect to the cutoff A. The shaded triangles repre¬ 
sent renormalized three-legged vertices with one bosonic and 
two fermionic external legs carrying the same sublattice label. 
For clarity, we have marked some of the external legs and the 
vertices with the associated sublattice labels. 



function e^{Q) = 1-1- The FRG flow of the 

polarization is given by 


'[CA(fc)+a(|fc-9l)P 


where we have now inserted the spin-degeneracy factor 
Ns = 2S + 1 = 2. It turns out that for K = Q = 0 the 
flow equation (5c) for the three-point vertices Tp{K, Q) 
reduces to the Ward identity [11] 7 aZa = 1, implying a 
partial cancellation between self-energy and vertex cor¬ 
rections in the above flow equations. Note that this 
cancellation is not properly taken into account in the 
random-phase approximation (RPA) where vertex cor¬ 
rections are assumed to be negligible [12]. Nevertheless, 
by combining the RPA with a RG procedure one can 
obtain accurate results for the renormalized velocity of 
graphene [13]. 

Although it is now straightforward to derive a closed 
system of FRG flow equations for Za and the two func¬ 
tions vt^{k) and let us neglect here the fre¬ 

quency dependence of the dielectric function, e{^{q,iuj) ~ 
eA{q). In this approximation Za = 7a = 1, but the renor¬ 
malization of the Fermi velocity is non-perturbatively 
taken into account. This is sufficient to obtain the correct 
quantum critical scaling in gaphene [10]. After perform¬ 
ing one of the integrations in Eqs. (9) and (10) we obtain 


AdAVA{k) 


A j"^ dip cos p 1 

‘2 k Jq tt y^l — 2(fc/A) cosp + (/c/A)2 eA(\/A2 — 2fcAcos(p -I- k^) ’ 


AdACAiq) = -2Nse^^ 



dp 0 ( 1 -I-^ cos</j — ^) sin^(/3 

TT ^[1 -b {q/2K)cospY - [9/2A]2 ['Ca(A) + {I + (q/A) cosp)va{A + qcosp)]' 


( ll a) 

( ll b) 


Note that Eq. (11b) has been obtained from Eq. (10) by 
shifting fe —> k + q/2 and then introducing elliptic coor¬ 
dinates. Eqs. (11a) and (lib) form a system of coupled 
integro-differential equations for the two momentum- and 
cutoff-dependent functions VAik) and CAiq)- The phys¬ 
ical renormalized velocity and the static dielectric func¬ 
tion are v{k) = limA->.o and e{q) = limA->.o eA( 9 )- 

We can easily recover the perturbative RG result (1) if we 
approximate ca ~ 1 on the right-hand side of Eq. (11a) 
and expand the integrand to leading order in k/A. How¬ 
ever, such an expansion is only valid for fc <C A, so that 
the physical limit A —>■ 0 at fixed A: 7 ^ 0 is not accessible 
within this approximation. 

We have solved the FRG flow equations (11a, lib) 
numerically without further approximations. Note that 
these equations are non-perturbative in the effective cou¬ 


pling constant a = jvp and should be quantitatively 
accurate even for large a. Our numerical result for VA{k) 
in Fig. 2 (a) clearly shows that the external momen¬ 
tum k replaces A as effective infrared cutoff as soon 
as A < fc. The physical momentum-dependent velocity 
v{k) = limA-j-o' 1 'a(A:) is shown in Fig. 2 (b). For compar¬ 
ison, we also show the velocities VA=k{k) and UA=fe( 0 ); 
the latter can also be obtained using the field-theoretical 
RG if one stops the RG flow at finite A and then substi¬ 
tutes A —>■ fc. While this recipe works perfectly to leading 
order in a where the perturbative result for u(fc) can be 
obtained by replacing A —)■ fc in Eq. (1), we see from Fig. 2 
(b) that this procedure remains approximately valid also 
for the physically relevant a = 2.2. Our FRG result for 
v(k) can be fitted by v(k)/vF = A(a) -I- B[a) ln(Ao/fc), 
where A(2.2) = 1.37 and B{2.2) = 0.51 for a = 2.2. 
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FIG. 2. (Color online) (a) Velocity VA{k) as a function of 
infrared cutoff A and external momentum k obtained from the 
numerical solution of Eqs. (11a) and (11b). In (b) the solid 
line represents our FRG result for the physical momentum- 
dependent velocity v{k) = VA=o{k) as a function of fc/Ao. For 
comparision, we also show VA=k(0) and VA=k{k). 


Note that our result for 5(2.2) is very close to the first 
order expression a/4 = 0.55. This is surprising, be¬ 
cause in a perturbative expansion the second order cor¬ 
rection leads to an additional contribution 620 ^ to the 
prefactor of ln(Ao/fc), where according to Mishchenko 
[16] &2 ~ —0.14, Vafek et al. [15] obtained 62 ~ —0.3, 
and Barnes et al. [14] found &2 ~ —0.32. In any case, for 
a = 2.2 the second order correction is substantial; our 
FRG calculation suggests that in this case the terms of 
order are to a large extent cancelled by higher order 
corrections. 

To compare our results with the experimental data for 
the renormalized velocity in suspended graphene [ 2 ], we 
need to fix the value of the ultraviolet cutoff Aq which 
appears as a free parameter in our low-energy action (3). 
Following the usual field-theoretical procedure [17], we 
eliminate Aq in favour of some measurable observable. 
Note that the authors of Ref. [2] use cyclotron resonance 
to measure the quasi-particle velocity v{kp) at the Fermi 
momentum kp for different densities n = kp/ir in slightly 
doped graphene. Assuming that the function v{k) is ap¬ 
proximately independent of n (which seems to be reason¬ 
able at low densities), we may fix Aq by demanding that 
our result for v{k = kp) agrees with the measured quasi¬ 
particle velocity at one particular value of fc^. Here we 
choose the data point at n « 49 x 10^°cm“^ to fix Ag; 
other choices lead to fits of similar quality. In Fig. 3 
we compare the renormalized velocity obtained with this 
prescription with the experimental data [2]. Obviously, 
in the entire range of available densities our FRG result 
agrees quite well with the data. 

Finally, in Fig. 4 (a) we show our numerical results 
for the momentum- and cutoff-dependent dielectric func¬ 
tion eA{q). Here the external momentum q and the cut¬ 
off A do not play the same role, because the dielectric 
function is defined in terms of the bosonic self-energy 
FIaIQ) while the infrared cutoff has been introduced in 
the fermionic propagator. The physical dielectric func- 


FIG. 3. (Color online) Comparison of our FRG result for 
the renormalized velocity v{kF) as a function of the density 
n = k% I TV with the data from Ref. [2] (dots with error-bars). 
The ultraviolet cutoff Aq is fixed as described in the text. 



FIG. 4. (Color online) (a) Momentum- and cutoff-dependent 
dielectric function ca)?) obtained from the numerical solu¬ 
tion of the FRG flow equations (11a) and (11b). (b) Physical 
dielectric function e{q) = limA-»o eA(</). The dashed line rep¬ 
resents the RPA result for Aq —>■ 00 . 


tion e{q) = €a=o{<i) is shown in Fig. 4 (b). Note that 
the logarithmic divergence of the velocity v{k) for small 
k leads to the logarithmic vanishing of the static bosonic 
self-energy n(( 7 ), so that the corresponding dielectric 
function e{q) logarithmically approaches unity for g —>■ 0 . 
Hence, in the static long-wavelength limit the Goulomb 
interaction in suspended graphene is not screened at all. 
On the other hand, if we replace the flowing velocity by 
the bare velocity on the right-hand side of our flow equa¬ 
tion (11b) and take the limit Ag — >■ 00 we recover the 
RPA result eRPA(<z) = 1 + 7 rA^sa/ 4 . 

In summary, we have derived and solved non- 
perturbative FRG flow equations for the momentum- and 
cutoff-dependent quasi-particle velocity UA(fc) and the 
static dielectric function eA(g) of suspended graphene. 
In the physical limit A —> 0 our result for v{k) diverges 
as ln(Ao/fc) for fc —0 and agrees very well with a recent 
experiment [2]. The dielectric function e{q) is shown to 
approach unity for g — 0 , in contrast to the prediction of 
the RPA. Our approach can be extended in several direc¬ 
tions: with some extra numerical effort, the frequency- 
dependence of the dielectric function can be taken into 
account. Moreover, we can generalize our approach to 
allow for the spontaneous formation of a charge-density 
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wave which transforms the system into an excitonic in¬ 
sulator [1, 18]. 

We thank A. Geim for making the data for the renor¬ 
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of us (C. B. and P. K.) acknowledge the hospitality 
of the physics department of the University of Florida, 
Gainesville, where part of this work has been carried out. 
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